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ABSTRACT 

Boutloukos et al. (2006) discovered twin-peak quasi-periodic oscillations (QPOs) in 1 1 observations of the 
peculiar Z-source Circinus X-l. Among several other conjunctions the authors briefly discussed the related 
estimate of the compact object mass following from the geodesic relativistic precession model for kHz QPOs. 
Neglecting the neutron star rotation they reported the inferred mass Mo = 2.2 + O.3M . We present a more 
detailed analysis of the estimate which involves the frame-dragging effects associated with rotating spacetimes. 
For a free mass we find acceptable fits of the model to data for (any) small dimensionless compact object 
angular momentum j = cJ/GM 2 . Moreover, quality of the fit tends to increase very gently with rising j. Good 
fits are reached when M ~ Mq[\ + 0.55(y + j 2 )]. It is therefore impossible to estimate the mass without the 
independent knowledge of the angular momentum and vice versa. Considering j up to 0.3 the range of the 
feasible values of mass extends up to 3M . We suggest that similar increase of estimated mass due to rotational 
effects can be relevant for several other sources. 

Subject headings: stars: neutron — X-rays: binaries 



1. INTRODUCTION 

Quasi-periodic oscillations (QPOs) appear in variabilities 
of several low-mass X-ray binaries (LMXBs) including those 
which contain a neutron star (NS). A certain kind of these os- 
cillations, the so-called kHz (or high-frequency) QPOs, come 
often in pairs with frequencies v L and v v typically in the range 
~50- 1300 Hz. This is of the same order as the range of fre- 
quencies characteristic for orbital motion close to a compact 
object. Accordingly, most kHz QPO models involve orbital 
motion in the inner regions of an accretion disk (see van der 
Kris 2006; Lamb & Boutloukos 2007, for a recent review). 

There is a large variety of QPO models related to NS 
sources (in some but not all cases they are applied to black 
hole (BH) sources too). Concrete models involve miscel- 
laneous mechanisms of producing the observed rapid vari- 
ability. One of the first possibilities proposed represents the 
"beat-frequency model assuming interactions between the 
accretion disk and spinning stellar surface (Alpar & Shaham 
1985; Lamb et al. 1985). Many other models primarily as- 
sume accretion disk oscillations. For instance, non-linear res- 
onance scenarios suggested by Abramowicz, Kluzniak and 
Collaborators (Abramowicz & Kluzniak 2001; Kluzniak & 
Abramowicz 2001; Abramowicz et al. 2003b, c; Horak 2008; 
Horak et al. 2009) are often debated. A set of the later mod- 
els join the beat frequency idea, magnetic field influence and 
presence of the sonic-point (Miller at al. 1998a; Psaltis et al. 
1999; Lamb & Coleman 2001, 2003). Some of numerous ver- 
sions of non-linear oscillation models and the late beat fre- 
quency models rather fade into the same concept that com- 
monly assumes the NS spin to be important for excitation of 
the resonant effects (Lamb & Coleman 2003; Kluzniak et al. 
2004; Petri 2005a,b,c; Miller 2006; Kluzniak 2008; Stuchlik 
et al. 2008; Mukhopadhyay 2009). Resonance, influence of 
the spin and magnetic field play a role also in the ideas dis- 
cussed by Titarchuk & Kent (2002) and Titarchuk (2002). An- 
other resonances are accommodated in models assuming de- 
formed disks (Kato 2007, 2008, 2009a,b; Meheut & Tagger 
2009). Further effects induced in the accreted plasma by the 
NS magnetic field (Alphen wave model, Zhang 2005; Zhang 



et al. 2007a,b), oscillations that arise due to comptonization 
of the disk-corona (Lee & Miller 1998) or oscillations ex- 
cited in toroidal disk (Rezzolla et al. 2003; Rezzolla 2004; 
Sramkova 2005; Schnittman & Rezzolla 2006; Blaes et al. 
2007; Sramkova et al. 2007; Straub & Sramkova 2009) are 
considered as well. At last but not least, already the kinemat- 
ics of the orbital motion itself provides space for consideration 
of "hot-spot-like" models identifying the observed variability 
with orbital frequencies. For instance recent works of Cadez 
et al. (2008) and Kostic et al. (2009) deal with tidal disruption 
of large accreted inhomogenities. Among the same class of 
(kinematic) models belongs also the often quoted "relativis- 
tic precession" (RP) kHz QPO model that is of our attention 
here. 

The RP model has been proposed in a series of papers by 
Stella & Vietri (1998, 1999, 2002). It explains the kHz QPOs 
as a direct manifestation of modes of relativistic epicyclic mo- 
tion of blobs arising at various radii r in the inner parts of the 
accretion disk. The model identifies the lower and upper kHz 
QPOs with the periastron precession v P and Keplerian v K fre- 
quency, 



v L (r) = v P (r) = v K (r) - v r (r), v v (r) = v K (r), 



(1) 



where v r is the radial epicyclic frequency of the Keplerian 
motion. Note that, on a formal side, for Schwarzschild space- 
time where v K equals a vertical epicyclic frequency this iden- 
tification merges with a model assuming m = — 1 radial and 
m = -2 vertical disk-oscillation modes. 

In the past years, the RP model has been considered among 
the candidates for explaining the twin-peak QPOs in several 
LMXBs and related constraints on the sources have been dis- 
cussed (see, e.g., Karas 1999; Zhang et al. 2006; Belloni et al. 
2007a; Lamb & Boutloukos 2007; Barret & Boutelier 2008a; 
Yan et al. 2009). While some of the early works discuss these 
constraints in terms of both NS mass and spin and include also 
the NS oblateness (Morsink & Stella 1999; Stella et al. 1999), 
most of the published implications for individual sources fo- 
cus on the NS mass and neglect its rotation. 

Two simultaneous kHz QPOs with centroid frequencies 
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of up to 225 (500) Hz have recently been also found 
by Boutloukos et al. (2006) in 11 different epochs of 
RXTE / Proportional Counter Array observations of the pecu- 
liar Z-source Circinus X-l. Considering the RP model they 
reported the implied NS mass to be M ~ 2.2M . The esti- 
mate was obtained assuming the non-rotating Schwarzschild 
spacetime and was based on fitting the observed correlation 
between the upper QPO frequency and the frequency differ- 
ence Av = v v - v L . In this paper, we improve the analysis of 
mass estimate carried out by Boutloukos et al. In particular, 
we consider rotating spacetimes that comprehend the effects 
of frame-dragging and fit directly the correlation between the 
twin QPO frequencies. We show that good fits can be reached 
for the mass-angular-momentum relation rather than for the 
preferred combination of mass and spin. 

2. DETERMINATION OF MASS 

Spacetimes around rotating NSs can be with a high pre- 
cision approximated by the three parametric Hartle-Thorne 
(HT) solution of Einstein field equations (Hartle & Thorne, 
1968; see Berti et al., 2005). The solution considers mass M, 
angular momentum J and quadrupole moment Q (supposed 
to reflect the rotationally induced oblateness of the star). It 
is known that in most situations modeled with the present 
NS equations of state (EoS) the NS external geometry is very 
different from the Kerr geometry (representing the "limit" of 
HT geometry for q = QM/J 2 — » 1). However, the situation 
changes when the NS mass approaches maximum for a given 
EoS. For high masses the quadrupole moment does not induce 
large differences from the Kerr geometry since q takes values 
close to unity (Appendix A.l). 

The previous application of the RP model mostly implied 
rather very large masses (e.g. Belloni et al. 2007a). These 
large masses are only marginally allowed by standard EoS. 
Also the mass inferred by Boutloukos et al. (2006) takes val- 
ues above 2M Q . Motivated by this we use the limit of two- 
parametric Kerr geometry to estimate the influence of the spin 
of the central star in Circinus X-l (see Appendix A.l where 
we pay a more detailed attention to rationalization and discus- 
sion of this choice allowing usage of simple and elegant Kerr 
formulae. 

2.1. Frequency Relations 

Assuming a compact object of mass Mcgs = GM/c 2 and 
dimensionless angular momentum j = cJ/GM 2 described by 
the Kerr geometry, the explicit formulae for angular veloci- 
ties related to Keplerian and radial frequencies are given by 
the following relations (see Aliev & Galtsov 1981; Kato et al. 
1998, or Torok & Stuchlfk, 2005) 



where T = c 3 /(2jtGM) is the "relativistic factor" and x = 
r/McGS- Considering Equations (1) and (2), we can write for 
v L and v v , both expressed in Hertz (see also Appendix A. 1.2 
where we discuss a linear expansion of this formula), 



T '- J v u ~\T-jv v 



2/3 



-3/ 



r-jvu 



4/3 



1/2 



(3) 



In the Schwarzschild geometry, where j = 0, Equation (3) 
simplifies to 



v L = V v < 1 



1-6 



1/2 



(4) 



leading to the relation 

Av = v v y]l-6(27TGMv LI ) 2/3 lc 2 (5) 
that was used by Boutloukos et al. for the mass determination. 

2.2. "Ambiguity" in M 

There is a unique curve given by Equation (3) for each dif- 
ferent combination of M and j (see Appendix A. 2 for the 
proof). The frequencies v L and v v scale as 1 /M and, as illus- 
trated in the left panel of Figure 1, they increase with growing 
j. Naturally, one may ask an interesting question whether for 
different values of M and j there exist some curves that are 
similar to each other. We investigate and quantify this task in 
Appendix A. 2. 

There we infer 1 that for j up to ~ 0.3 one gets a set of nearly 
identical integral curves where M, j and Mo roughly relate as 
follows 



with 



M=[l+k(j + f)]M 



= 0.7 . 



(6) 



This result is illustrated in the right panel of Figure 1 . Clearly, 
when using relation (6), any curve plotted for a rotating star 
of a certain mass can be well approximated by those plotted 
for a non-rotating star with a smaller mass, and vice versa. 
Furthermore we find that (see Appendix A. 2) when the top 
parts of the curves (corresponding to v v /v L ~ 1-1.5 ) are 
considered only, the best similarity is reached for 

k = 0.75 . 

These parts of the curves are potentially relevant to most of the 
atoll and high-frequency Z-sources data. On the other hand, 
for the (bottom) parts of the curves that are potentially rele- 
vant to low-frequency Z-sources including Circinus X-l, the 
best similarity is achieved for 

k = 0.65 (0.55, 0.5) when v L ,/v L ~ 2 (3, 4) . 

Taking into account the above consideration we can expect 
that the single-parameter best fit to the data by relation (4) 
roughly determines a set of mass-angular-momentum combi- 
nations (6) with similar x 2 . The result of Boutloukos et al. 
then implies that good fits to their data, displaying v v jv L ~ 3, 
should be reached for M ~ 2.2M [1 + 0.55(; + /)]. In what 
follows we fit the data and check this expectation. 



1 We first consider a special set of apparently similar curves sharing the 
terminal points. The set is (numerically) given by the particular choice of M, 
for any j implying the same orbital frequency at the marginally stable circular 
orbit. The curves then only slightly differ in their concavity that increases 
with growing j. 
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Figure 1. — Left: relation between the upper and lower QPO frequency following from the RP model for the mass M = 2.5M Q . The consecutive curves differ 
in j e (0, 0.3) by 0.05. Right: relations predicted by the RP model vs. data of several NS sources. The curves are plotted for various combinations of M and j 
given by Equation (6) with k = 0.7. The datapoints belong to Circinus X-l (red/yellow color), 4U 1636-53 (purple color) and most of other Z- and atoll- sources 
(black color) exhibiting large population of twin-peak QPOs. 
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Figure 2. — Left: x dependence on the parameters M and j assuming Ken' solution of Einstein field equations. The continuous white curve indicates the 
mass-angular-momentum relation (7). The continuous thin green curve denotes j giving the best^ 2 for a fixed M. The dashed and thick green curve indicates 
the same dependence but calculated using formulae (A2) and (A6) linear in /, respectively. The reasons restricting the calculation of the thick curve up to j = 0.4 
are discussed in Section A. 1 .2. Right: related profile of the best x 2 for a fixed M. The arrow indicates increasing j. 



2.3. Data Ma tch ing 

In the right panel of Figure 1, we show the twin-peak fre- 
quencies measured in the several atoll and Z-sources together 
with the observations of Circinus X-l. For the Circinus X- 
1 data we search for the best fit of the one-parametric rela- 
tion (4). Already from Figure 1, where these data are em- 
phasized by the red/yellow points, one may estimate that the 
best fit should arise for Mo e 2-2.5M . Using the stan- 
dard least squares method (Press et al. 2007) we find the 
lowest^- 2 = 15 = 2dof for the mass Mq = 2.2M which 
is consistent with the value reported by Boutloukos et al. 
The symmetrized error corresponding to the unit variation 
of x 1 ls ±O.3M . The asymmetric evaluation of Mo reads 
2.2[+0.3; -O.1]M . The white curve in Figure 2 indicates the 
mass-angular-momentum relation implied by Equation (6), 



M = 2.2M [1 +k(j + f)], £ = 0.55. 



(7) 



For the exact fits in Kerr spacetime we calculate the rel- 
evant frequency relations for the range of M € 1-4M and 
j e0-0.5. These relations are compared to the data in order to 
calculate the map of x 2 - We use the step equivalent to thou- 

2 After Barret et al. (2005a,b); Boirin et al. (2000); Belloni et al. (2007a); 
di Salvo et al. (2003); Homan et al. (2002); Jonker et al. (2002a,b); Mendez 
& van der Klis (2000); Mendez et al. (2001); van Straaten et al. (2000, 2002); 
Zhang etal. (1998). 



sand points in both parameters and obtain a two-dimensional 
map of 10 6 points. This color-coded map is included in the 
left panel of Figure 2. One can see in the map that the ac- 
ceptable x 2 is rather broadly distributed. The thin solid green 
curve indicates j corresponding to the best x 2 for a fixed M. 
It well agrees with the expected relation (7) denoted by the 
white curve. The right panel of Figure 2 then shows in de- 
tail the dependence of the best x 2 for the fixed M. It is well 
visible that the quality of the fit tends to very gently, mono- 
tonically increase with rising j and it is roughly^- 2 ~ 15 for 
any considered j. 

3. DISCUSSION AND CONCLUSIONS 

The quality of the fit tends to very gently, monotonically 
increase with rising j and it is roughly 

X 2 ~2dof o M ~ 2.2[+O.3,-O.l]M x[l+O.55(;+; 2 )]. (8) 

Therefore one cannot estimate the mass without the indepen- 
dent knowledge of the spin or vice versa, and the above re- 
lation provides the only related information implied by the 
geodesic RP model. 

To obtain relation (8), the exact Kerr solution of Einstein 
field equations was considered. The choice of this two- 
parametric spacetime description and related formulae (2) is 
justified by a large value of the expected mass Mo (see Ap- 
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pendix A.l for details). In Appendix A. 1.2 we discuss the 
utilization of the linearized frame-dragging description. Fig- 
ure 2 includes the mass-spin dependence giving best x 1 re- 
sulting when the fitting of datapoints is based on the associ- 
ated formulae (A2) respectively (A6). Considering that v L (vu) 
formula (3) merge up to the first order in j with the v L (v v ) re- 
lation (A6) linear in j one can expect that the associated M{ j) 
relations obtained from fitting of data should roughly coin- 
cide up to j ~ 0.1-0.2. From the figure we can find that 
there is not a big difference between the resulting M(j) rela- 
tions even up to much higher j. The extended coincidence can 
be clearly explained in terms of the kHz QPO frequency ratio 
R = Vu/v L . 3 

Observations of Circinus X-l result to R ~ 2.5-4.5 while 
usually it is R ~ 1.2-3 (and most often R ~ 1.5, Abramow- 
icz et al. 2003b; Torok et al. 2008b; Yan et al. 2009). As- 
suming the RP model along with any j e (0, 1), the ratio 
R — 2 corresponds with a good accuracy to radii where the 
radial epicyclic frequency reaches its maximum (Torok et al. 
2008a). Only values lower than R ~ 2 are then associated 
with the proximity of innermost stable circular orbit (ISCO) 
where the effects of frame dragging become to be highly non- 
linear in both j and r. Accordingly, for a given j, in the case 
when R ~ 3, the individual formulae restricted up to certain 
orders in j are already close to their common linear expansion 
in j and differ much less than for R ~ 1 .5 (see Appendix A. 1). 

The rarely large R and associated high radial distance (both 
already remarked by Boutloukos et al. 2006, although in a 
different context) in addition to large Mq warrant relevance 
of relation (8) for rather high values of the angular momen- 
tum. Consequently, we can firmly conclude that the upper 
constrained limit of the mass changes from the value 2.5M 
to 3M for j = 0.3 and even to 3.5M for j = 0.5. The Value 
of Mo that is above 2M and the increase of M with grow- 
ing j for corotating orbits elaborated here are challenging for 
the adopted physical model. Further detailed investigation in- 
volving realistic calculations of the NS structure can be there- 
fore effective in relation to EoS selection or even falsifying 
the RP model. 

Finally, we note that discussed trend of increase of esti- 
mated mass arising due to rotational effects should be rel- 
evant also for several other sources. Of course, many sys- 
tems display mostly low values of R. These low values of 
R are in context of the RP model suggestive of proximity of 
ISCO. Torok (2009) and Zhang et al. (2010) pointed that un- 
der the consideration of the RP model and j = 0, most of the 
high-frequency sources data are associated with radii close to 
r = 6.75M. Possible signature of ISCO in high frequency 
sources data has been also reported in a series of works by 
Barret et al. (2005a; 2005b; 2006) based on a sharp drop in 
the frequency behavior of the kHz QPO quality factors (for 
instance the atoll source 4U 1636-53 denoted by "blueberry" 
points in Figure 1) clearly exhibits both low R and a drop of 
QPO coherence, see Boutelier et al. 2009). Considering the 
proximity of ISCO, high-order non-linearities in both j and 
r are important and even small differences between the actual 
NS and Kerr metric could have certain relevance. For this rea- 
son some caution is needed when applying our results to high 
frequency sources. 

3 Orbital frequencies scale with l/M. For any model considering v L and 
v u given by their certain combination, the ratio R represents the measure of 
radial position of the QPO excitation (provided that the NS spin and EoS are 
fixed). 
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APPENDIX 

A. APPROXIMATIONS, FORMULAE AND 
EXPECTATIONS 

A.l. Matching Influence of Neutron Star Spin 

Rotation and the related frame-dragging effects strongly in- 
fluence the processes in the vicinity of compact objects and 
there is a need of their reflection in the appropriate spacetime 
description. External metric coefficients related to up-to-date 
sophisticated models of rotating NS are taken out of the model 
in two distinct ways. In the first way, the coefficients are ob- 
tained "directly" from differential equations solved inside the 
numeric NS model, while in the second (more usual) way 
they are inferred from the main parameters of the numeric 
model (mass, angular momentum etc.) through an approxi- 
mative analytic prescription. Several commonly used numeri- 
cal codes related to rotating NS have been developed and dis- 
cussed (see, RNS: Stergioulas & Morsink, 1997; LORENE: 
Gourgoulhon et al., 2000; and also Nozawa et al., 1998; Ster- 
gioulas & Friedman, 1995; Cook et al., 1994; Komatsu et al., 
1989). 

A. 1.1. Analytical Approximations and High-mass Neutron Stars 

In the context of a simplified analysis of NS frame-dragging 
consequences, an approximation through two solutions of 
Einstein field equations is usually recalled: Lense-Thirring 
metric also named linear-Hartle metric (Thirring & Lense 
1918; Hartle & Sharp 1967; Hartle 1967) and Kerr-black-hole 
metric together with related formulae (Kerr 1963; Boyer & 
Lindquist 1967; Carter 1971; Bardeen et al. 1972). It is ex- 
pected that the Lense-Thirring metric well fits the most im- 
portant changes (compared to the static case) in the external 
spacetime structure of a slowly rotating NS. This expectation 
is usually assumed for / < 0.1 -0.2. 4 Due to asymptotical 
flatness constraints the formulae related to Lense-Thirring, 
Kerr and some other solutions considered for rotating NS 
merge when truncated to the first order in j. Accordingly, 
for the astrophysical purposes there is a widespread usage of 
the approximative terms derived with the accuracy of the first 
order in j. While these approximations are two-parametric, 
the more realistic approximations; for instance those given 
by the HT metric (Hartle & Thorne 1968) and related terms 
(Abramowicz et al. 2003a), relations of Shibata & Sasaki 
(1998) or the solution of Pachdn et al. (2006), deal with more 
parameters and provide less-straightforward formulae. Per- 
haps also because of that they are not often considered in dis- 
cussions of concrete astrophysical compact objects. 

4 The interval < j < 1 X 10~' is often assumed as one of the several 
possible definitions of "slow rotation". However, in relation to implications 
of the frame-dragging effects, the effective size of this interval depends on the 
radial coordinate. For x close or below ,v ras the interval in j rather reduces to 
low values. On the other hand for x above the radius of the maximum of v r 
the interval can be extended to j higher than j = 0.2. The term slow rotation 
is also frequently considered in another context. For instance when using 
the HT metric in NS models the slow rotation is usually associated with the 
applicability of the metric and consequently to spins up to ~ 800Hz for most 
of EOS and NS masses. For these reasons we do not use the term elsewhere 
in the paper. 




Figure 3. — Left: parameter q for several EoS. Shaded areas denote q = 6 and q = 3. — Right: ISCO frequencies for the same EoS as used in the left panel. 
The curves are calculated for mass 1 AM e and a relevant maximal allowed mass. The curves following from the exact Kerr solution and linear relation ( A4) are 
displayed as well. The quadratic relation denoted by the black-dashed curve is discussed later in Section A.2. 1 . 
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x x 

Figure 4. — Frequencies of the perturbed circular geodesic motion. Relations for the Kerr metric given by Equation (A2) are denoted by blue and dashed-blue 
curves. Relations (A2) are indicated by red curves while relation (A5) is plotted using the green color. Dotted relations denote the Kerr- and linearized- vertical 
frequencies that are not discussed here (see Morsink & Stella 1999; Stella et al. 1999). Inset emphasizes a difference between the radii fulfilling the ISCO 
condition v r = for the relations (2) explicitly given by Equation (Al), Equation (A2) and the ISCO-radius given by Equation (A3). 




(M/M Q ) x v v [kHz] j 

Figure 5. — Left: the RP model frequency relations given by Equation (3), blue curves; formulae (A2), red curves; relations (A6), green curves. Relation (A7) 
roughly determining the applicability of (A6) is denoted by the dashed black/yellow curve. Right: related differences Av between the lower QPO frequency 
implied by the Kerr formulae (3) and those following from Equation (A2) respectively Equation (A6) indicated by continuous respectively dashed curves. 
Different colors correspond to different frequency ratio R. Shaded areas indicate Av < 5% and Av < 2%. 



Astrophysical applicability of the above analytical ap- 
proaches has been extensively tested in the past ten years. Cri- 
teria based on the comparison of miscellaneous useful quan- 
tities have been considered for these tests (e.g., Miller et al. 
1998b; Berti et al. 2005). It has been found that spacetimes 



induced by most of up-to-date NS EoS without inclusion of 
magnetic field effects are well approximated with the HT so- 
lution of the Einstein field equations (see Berti et al. 2005, for 
details). The solution reflects three parameters: NS mass M, 
angular momentum J and quadrupole moment Q. Note that 



6 



TOROK ET AL. 



Kerr geometry represents the "limit" of the HT geometry for 
q = Ql J 2 — > 1. The parameter q then can be used to charac- 
terize the diversity between the NS and Kerr metric. 

The left panel of Figure 3 displays a dependence of q on 
the NS mass. This illustrative figure was calculated follow- 
ing Hartle (1967), Hartle & Thorne (1968), Chandrasekhar & 
Miller (1974) and Miller (1977). The considered EoS are de- 
noted as follows (see Lattimer & Prakash, 2001 and Lattimer 
& Prakash, 2007 for details): 
[EoSl] SLy 4, Rikovska Stone et al. (2003). 
[EoS2] APR, Akmal et al. (1998). 

[EoS3] AU (WFF1), Wiringa et al. (1988); Stergioulas & 
Friedman (1995). 

[EoS4] UU (WFF2), Wiringa et al. (1988); Stergioulas & 
Friedman (1995). 

[EoS5] WS (WFF3), Wiringa et al. (1988); Stergioulas & 
Friedman (1995). 

Inspecting the left panel of Figure 3 we can see that for EoS 
configurations resulting to low- or middle-mass of the central 
star (M up to 0.8M max , i.e., roughly up to 1.4M- 1.8M ) de- 
pending on EOS, the implied HT geometry is rather different 
from the Kerr geometry. More specifically, for a fixed central 
density, q strongly depends on the given EoS and substantially 
differs from unity. On the contrary, for high mass configura- 
tions q approaches unity implying that the actual NS geometry 
is close to Kerr geometry. One can expect that in such case 
formulae related to the Kerr geometry should provide better 
approximation than for low values of M. Next, focusing on 
high-mass NS, we briefly elaborate some points connected to 
the applicability of the Kerr formulae and related linearized 
terms. 

A. 1.2. Kerr and linearized Kerr Formulae: Comparison, 
Utilization and Restrictions 

The radial epicyclic frequency goes to zero on a particular, 
so-called marginally stable circular orbit x ms (e.g., Bardeen 
et al. 1972). In Kerr spacetimes it is given by the relation 
(Bardeen et al. 1972) 



x ms = 3 + Z 2 - V(3 - Zi)(3 + Zi + 2Z 2 ) , (Al) 
where 

Zl = i + (i-f)^[(i +J f' 3 + (i-jy' 3 ], 



+ Z 2 



Below x m there is no circular geodesic motion stable with re- 
spect to radial perturbations. The orbit is often named ISCO 
and determines the inner edge of a thin accretion disk. The 
corresponding ISCO orbital frequency v s (x ms ) represents the 
highest possible orbital frequency of the thin disk and the re- 
lated "spiraling" inhomogenities (Kluzniak et al. 1990). De- 
pendence of ISCO frequency on j following from Equation 
(Al) is shown in the right panel of Figure 3. 

Assuming the description of geodesic motion accurate in 
the first order of j, using Taylor expansion around j = 0, one 
may rewrite the explicit terms in Equation (2) as 



ni i - - + 



8./ 



(A2) 



Consequently, linearized formula for the ISCO radius can be 
expressed as 



(A3) 



*ms = 6 - 4 A / - j . 



Note that the root of the expression for to 2 from Equation (A2) 
is of higher order in j so that the exact radius where w r van- 
ishes agrees with the solution (A3) only in the first order of j. 
The related ISCO frequency can be evaluated as (Kluzniak & 
Wagoner 1985; Kluzniak et al. 1990) 

v K (x ms ) = (MJM) x (1 + 0.749;") x v,(x ms , M = M e , j = 0) 
= (M Q /M) x (1 + 0.749 j) x 2197 Hz. (A4) 

This frequently considered relation is included in the right 
panel of Figure 3. 

In the right panel of Figure 3 we integrate the ISCO fre- 
quencies plotted for several EoS (the same as in the left panel). 
We choose two groups of models - one calculated for the set 
of five different EoS and "canonic" mass 1 AM Q , the other one 
for the same set of EoS but considering a maximal mass al- 
lowed by each individual EoS. This choice allows for the illus- 
tration of middle- and high- mass behavior of ISCO relations 
and comparison of their simple approximations. Clearly, for 
middle mass configurations Equation (A4) provides better ap- 
proximation than using the Kerr-spacetime formulae (see also 
Miller et al. 1998b). On the other hand, when the high-mass 
configurations are considered, the Kerr solution provides bet- 
ter approximation than Equation (A4). Moreover, its accu- 
racy is higher than the accuracy of both approximations for 
middle-mass configurations. 

A. 1.3. Geodesic frequencies and RP model 

When the expression for the radial epicyclic frequency 
given by Equation (2) or (A2) is fully linearized in j, it leads 
to 

(jc - 6)x 3 ' 2 + 3 j(x + 2) 



VC* - 6)x 7 



(A5) 



Relation (A5) provides a good approximation except for the 
vicinity of x ms (j) as it diverges at x = 6. Note that this diver- 
gence arises only for corotating but not counterrotating orbits 
(which we however do not discuss in this paper). For any pos- 
itive j < 0.5 the fully linearized frequency (Equation (A5)) 
does not differ from u> r given by Equation (A2) for more than 
about 5% when x > 6 + Aj. The left panel of Figure 4 com- 
pares the frequencies of geodesic motion associated directly 
with Kerr metric, to formulae (A2) and (A6), respectively. 

Assuming linearized Keplerian frequency given by Equa- 
tion (A2) and the radial epicyclic frequency (Equation (A5)), 
we can write for the RP model the relation between v L and v v 
as 



1 - Vl -6a + 



2jv„ (a - 2) 



IFV1 -6a 



(A6) 

which equals the first-order expansion of Kerr spacetime 
Equation (3) and also to the first-order expansion of the same 
relation if it would be derived for Lense-Thirring or HT met- 
ric. Similarly to relation (A5), relation (A6) loses its physical 
meaning for frequencies close to v K (ISCO) since it reaches a 
maximum at frequencies that can be expressed with a small 
inaccuracy as 



12-v„/200\ M 



)(Hz). 



(A7) 



The left panel of Figure 5 compares the frequency relations 
(A6) to relations (3) and those following from formulae (A2). 
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It is useful to discuss their differences in terms of the fre- 
quency ratio R = v v /v L . For a fixed j the frequencies v L and 
v v scale with 1 /M. The ratio R then represents a "measure" 
of the radial position of the QPO excitation. It always reaches 
R = 1 at ISCO where the non-linear j terms are important and 
R — at infinity where the spacetime is flat. Note that R — 2 
almost exactly corresponds to the maximum of v r for any j 
(Torok et al. 2008a). 

The right panel of Figure 5 quantifies differences between 
the QPO frequency implied by the Kerr formulae (3), relations 
(A2) and relation (A6). We can see that differences between 
the Kerr relations (3) and those implied by formulae (A2) be- 
come small when R > 2 (Av < 5% for j < 0.5). For/? ~ 3 and 
higher, relations (3) and those implied by Equation (A2) are 
almost equivalent nearly merging to their common linear ex- 
pansion (A5). Note that taking into account relation (A7) the 
linear expansion (Equation (A6)) provides reasonable physi- 
cal approximation for spins and frequency ratios roughly re- 
lated as 

j< O.30R-1). (A8) 



A. 1.4. Applications 

Several values of NS mass previously reported to be re- 
quired by the RP model, including the estimate of Boutloukos 
et al., belong to the upper part of the interval allowed by stan- 
dard EoS. We can therefore expect low q and take the advan- 
tage of the exact Kerr solution for most of the practical cal- 
culations needed through the paper. Unlike to formulae trun- 
cated to certain order, all the formulae derived from the exact 
Kerr solution are from the mathematical point of view fully 
selfconsistent for any j. This allows us to present the content 
of Appendix A. 2 in a compact and demonstrable form. 

In Section 3 we finally compare the results of QPO fre- 
quency relation fits for Circinus X-l using the Kerr solution 
and those done assuming Equations (A2) respectively (A6). 
From the previous discussion it can be expected that for Circi- 
nus X-l, due to its exceptionally high R, the fits obtained 
with the Kerr formulae (3) and "linear" formulae (A2) should 
nearly merge with the fits obtained assuming the common lin- 
ear expansion (Equation (A6)). Note also that, on a technical 
side, the linear expansion can be used up to j ~ 0.3 - 0.4 since 
the lowest R in the Circinus X-l data is R ~ 2-2.5 (Equation 
(A8)). 

A.2. Uniqueness of Predicted Curves and "Ambiguity" in 

M 

The radial epicyclic frequency vanishes at x ms . In the RP 
model it is then vj™* = v t max = v A (jic ms ). Obviously, if there 
are two different combinations of M and j which, based on 
the RP model, imply the same curve v v (v L ), such combinations 
must also imply the same ISCO frequency. 

In the left panel of Figure 6 we show a set of curves con- 
structed as follows. We choose M* = 2.5M and j e (0, 0.5) 
and for each different j we numerically find M such that the 
corresponding ISCO frequency is equal to those for M* Q and 
j = 0. Then we plot the v v {v L ) curve for each combination of 
M and j. We can see that except for the terminal points the 
curves split. The frequencies in the figure can be rescaled for 
any "Schwarzschild" mass Mo as M^/Mq. Thus, the scatter 
between the curves provides the proof that one cannot obtain 
the same curve for two different combinations of M and j. 



On the other hand, the discussed scatter is apparently small 
and the curves differ only slightly in the concavity that grows 
with increasing j. This has an important consequence. The 
curves are very similar with respect to the typical inaccuracy 
of the measured NS twin-peak data and there arises a possible 
mass-angular-momentum ambiguity in the process of fitting 
the datapoints. Next, we derive a simple relation approximat- 
ing this ambiguity. 

A.2.1. Formulae for ISCO Frequency 

The ambiguity recognized in the previous section is implic- 
itly given by dependence of the ISCO frequency on the NS 
angular momentum which for the Kerr metric follows from 
relations (2) and (Al). In principle we can try to describe 
the ambiguity starting with these exact relations. The other 
option is to assume an approximative formula for the ISCO 
frequency. One can expect that this formula should be at least 
of the second order in j if consideration of spin up to j = 0.5 
is required. We check an arbitrarily simple form 

v A .(x ms ) = (M /M) x[l+k(j + /)] x 2197 Hz. (A9) 

The right panel of Figure 6 indicates the square of difference 
between the exact ISCO frequency in Kerr spacetimes follow- 
ing from Equation (Al) and the value following from Equa- 
tion (A9). Inspecting the figure we can find that the particular 
choice of k — 0.75 provides a very good approximation. 

Figure 7 then directly compares the exact relation and re- 
lation (A9) with k = 0.75. For comparison, the first-order 
Taylor expansion formula (A4) is indicated. Clearly, using 
Equation (A9) one may well approximate the Kerr-ISCO fre- 
quency up to j ~ 0.4 and describe the discussed ambiguity in 
terms of Schwarzschild mass Mq as 

M~[l+k(j + f)]M , (A10) 

where k = 0.75. In further discussion we therefore assume 
this formula. 

A.2.2. Comparison Between Curves 

The curves given by Equation (A10) with k = 0.75 are il- 
lustrated in the left panel of Figure 8. Here we quantify their 
(apparent) conformity and investigate its dependence on k. It 
is natural to consider the integrated area S between the curve 
for Mo, j ' - and the others as the relevant measure. The 
right panel of Figure 8 shows this area as the function of k 
in Equation (A10) for several values of j. The same panel 
also indicates the values related to the set of curves for mass 
found numerically from the exact Equations (Al), i.e., curves 
in the left panel of Figure 6. We can see that values of S for 
k = 0.75 are comparable to those related to Figure 6. More- 
over, for a slightly different choice of k — 0.7, all the values 
are smaller. The ambiguity in mass with relation (A10) is 
therefore best described for k ~ 0.7 when the data uniformly 
cover the whole predicted curves. 

The available data are restricted to certain frequency ranges 
and often exhibit clustering around some frequency ratios 
v L ,/v L (see Abramowicz et al. 2003b; Belloni et al. 2007b; 
Torok et al. 2008a; Torok et al. 2008b,c; Barret & Boute- 
lier 2008b; Torok 2009; Boutelier et al. 2009; Bhattacharyya 
2009). It is then useful to separately examine the mass ambi- 
guity for related segments of the curves. Such investigation is 
straightforward for small segments. Let us focus on a single 
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Figure 7. — Left: ISCO frequency calculated from Equation (A9) vs. exact relation implied by the Kerr solution (dashed vs. thick curve). The linear relation 
(A4) is shown as well for comparison (dotted curve). Right: the related relative difference from ISCO frequency in Kerr spacetime. 
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Figure 8. — Left: the set of curves plotted for combinations of M and j given by Equation (A 10) with k = 0.75. Right: the integrated area S related to Equation 
(A10). Different values of j are color-coded. The same color-code is relevant for horizontal lines. These lines denote the values of S arising for the set of curves 
numerically found from Equation (A 1 ) and plotted in the left panel of Figure 6. The two red vertical lines denote the case of k = 0.75 (curves v u (v L ) shown in the 
left panel of this figure) respectively k = 0.7 (see the text for explanation). 




point [v L , Vu] representing a certain frequency ratio for a non- 
rotating star ( j = 0) of mass Mq. Assuming relation (A10) one 
may easily calculate the value of k which rescales the mass to 
M + Mq for a fixed non-zero j in order to get exactly the same 
point [v L , v v ]. We applied this calculation for v v /v L e (1, 10) 
and j € (0, 0.5). The output is shown in Figure 9. From 
the figure, it is possible to find k that should best describe the 
ambiguity for a given frequency ratio (and thus for a small 



segment of data close to the ratio). It also indicates the length 
of the curve v v (v L ) integrated from the terminal (i.e., ISCO) 
point to the relevant frequency ratio (assuming j = 0). This 
length is given in terms of the percentual share / on the total 
length £ of the curve Vu(v L ), whereas the absolute numbers 
scale with 1 /M. 

Apparently, the segments with v v /v L € (1, 2) cover about 
50% of the total length £. while k only slightly differs from 
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TABLE 1 

The coefficient k representing mass-angular momentum ambiguity (A10). 

segment k in M ~ [1 + k(j + j z )]M I [%] distance from ISCO xM /M [km] 



vjvy ~ 1.5 0/75 23 T 

v L /v u ~ 2 0.65 50 3 

vjv v ~ 3 0.55 70 7 

Vi/v K ~ 4 0.50 80 12 

vjvu ~ 5 0.45 83 16 

yjyj, ~ 6 O40 85 20 

whole curve 077 



the value of 0.7. We recall that this top part corresponds to 
the most of atoll and Z-sources data. For the segments of 
curves related to the sources exhibiting high frequency ratios 
such as Circinus X-l, there is an increasing deviation from 
the 0.7 value and the coefficient reaches k ~ 0.6-0.5. A more 
detailed information is listed up to v v /v L = 5 in Table 1 pro- 
viding the summary of this section. 



I [%] 

50 70 80 85 90 95 

1.0 [ 1 



0.8 : 




o r 

2 4 6 8 10 

Vu/V L 

Figure 9. — Values of k approximating the M — j ambiguity for the individ- 
ual segments. The upper axes indicate the length of the curve v u (y L ) integrated 
from the ISCO point to the relevant frequency ratio. 
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